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I. INTRODUCTION 



One of the central issues in the field of high temperature superconductivity is the symmetry of the superconducting 
order parameter. A number of recent experiments, specifically Josephson interference experiments, suggest a d-wave 
order parameter with d x i_ y i symmetry,Efu while experiments on c-axis tunnelinga are consistent with a traditional 
s-wave order parameter. In this paper we calculate the in-plane and c-axis conductivities for layered, cZ-wave supercon- 
ductors in the limit q — > and for frequencies fiuj ~ A (infrared region). This long wavelength regime is appropriate 
for superconductors with a mean free path that is short compared to the penetration length, I <C A. Given the small 
coherence length in cuprates, the q — ► limit for the conductivity easily accommodates the clean limit condition, 
| £ <C I, which is necessary for d-wave superconductivity not to be suppressed. Calculations of the conductivity in the 
ultraclean limit > A > (), which include polarization sensitivity to the in-plane absorption and vertex corrections 
' to the conductivity from collective modes, will be discussed in a separate paper. 
ON ■ In order to represent the layered structure of high-T c superconductors, we consider the microscopic model of 
coupled superconducting layers introduced in Ref. [| In this model the charge carriers of each layer form a two- 
dimensional Fermi liquid, which we describe by the quasiclassical theory of superconductivity. Thus, in-plane transport 
is dominated by charged quasiparticles moving with an in-plane Fermi velocity, v/. Interiayer transport may originate 
from coherent prQpagation,B~a static ancLthcrmally activated (e.g., phonon assisted QLuair fluctuation) incoherent 
charge transfer resonant transfer Jiil or several transfer mechanisms in parallelOcj Interiayer transport has 
also been argued to take place through coherent pair tunneling.EU In this report we consider a model in which the 
interiayer coupling is mediated through incoherent hopping processes (the interiayer diffusion model), as in Ref. [(]. 
These incoherent scattering events are predominantly responsible for charge transfer in the c direction and Josephson 
coupling between the planes. For simplicity, we consider a circular Fermi surface for each conducting layer. This 
model is expected to be appropriate for superconductors whose c-axis behavior is that of a stack of superconductor- 
insulator-superconductor (SIS) Josephson junctions. Such behavior has recently been observed experimentally in 
several layered superconducting compounds .E3c3 

The infrared conductivity of three-dimensional bulk superconductors with unconventional pairing has been cal- 
culated previously by KJemm et al.jEl Hirschfeld et al.,t3 and for two-dimensional superconductors by Hirschfeld, 
Putikka and ScalapinccS using Green function methods; our results agree with those found previously whenever a 
direct comparison is possible. We calculate the infrared conductivity for layered superconductors and find striking 
differences between the conductivities of the d-wave and s-wave models. These differences between s-wave and d- 
wave conductivities are not due just to the differences in the density of states, which are known to lead to different 
low-temperature behavior of the penetration depth, NMR relaxation rates, etc.CJ There are additional features in 
the conductivity of d-wave superconductors which are attributed to the formation of optically active Andreev bound 
states. 

The outline of this paper is the following . In Sec. |J, we present our microscopic model in terms of the Fermi liquid 



theory of superconductivity. Then in Sec. Ill we formulate the linear current response to a spatially homogeneous, 
electromagnetic perturbation. Section IV includes a discussion of the phenomenological parameters of our model as 
well as the results of the current response calculations. We relate the model parameters to measurable normal-state 
quantities and calculate properties of the superconducting state for both s-wave and <i-wave pairing interactions. We 
discuss the consistency between these two models and experiments, and discuss the results of the linear response 
calculations of the in-plane and c-axis infrared conductivities. 
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II. QUASICLASSICAL THEORY 



We perform the calculations within the quasiclassical theory of superconductivity, first developed by Eilenberger. El 
Larkin and OvchinnikovJ22l and Eliashberg.E3 This theory is capable of describing both the equilibrium and dynamical 
properties of a Fermi liquid, in either the normal or superconducting state. We briefly outline the basic equations 
of this theory without discussing their justification^, readers interested in these details should see the original papers 
mentioned above, or the recent review articles .EjED We follow closely the notation of Refs. ^| and |33|. In this report 
we are interested in the linear response of a homogeneous, superconducting Fermi liquid to an electric field described 
by a spatially uniform, time-dependent vector potential. With this in mind, we start from the basic equations of the 
quasiclassical theory neglecting the dependence on the space coordinate (R). 



A. Basic Equations 



The most fundamental quantity in the quasiclassical theory of superconductivity is the matrix propagator for 
quasiparticle excitations. In-order to describe noncquilibrium properties we use Keldysh's formulation of dynamical 
phenomena (see Appendix)B3 This technique requires three types of propagators: retarded, g R , advanced, g A , and 
Keldysh, g . These propagators are 2x2 matrices in particle-hole (Nambu) space, 



g R < A > K (s;e,t)=( %* K i R , K I. (1) 



gR.,A,K fR,A,K 
jR,A,K gR,A,K 



Note that we are limiting our considerations to distribution functions which are spin_scalars only; a more complete 
description of this theory, including spin-dependent effects, can be found elsewhereo For a spatially homogeneous 
Fermi liquid, the propagators g R ' A > K (s; e, t) are functions of the Fermi surface position, s, the single particle excitation 
energy, e, and time, t. Their diagonal components are related to the spectrum and distribution of Bogoliubov 
quasiparticle excitations, while their off-diagonal components yield the Cooper-pair amplitude. In principle, g K 
contains all of the information about the relevant measurable quantities (equilibrium and dynamic) in both the 
normal and superconducting states. We make use of a compact notation which combines the three Nambu matrices, 
gR.A.K ^ j n ^ a sm gi e Nambu-Keldysh matrix, 

9(s;e,t)=( § * J*Y (2) 

Thus, any operator written without superscript (e.g., a) is interpreted as a Nambu-Keldysh matrix, while its retarded, 
advanced, and Keldysh Nambu-matrix components are denoted by the appropriate superscripts (e.g., a R ' A ' K ). 
The central equation of the Fermi liquid theory of superconductivity is the quasiclassical transport equation 



c 

supplemented by the normalization condition 



v/ ■ A e ) f 3 - e$ £ l - A e - ±i , gi 



0, (3) 



k®ge. = -vri. (4) 



The f 3 in equation (Jfy is to be understood as the third Pauli matrix in particle-hole space combined with the unit 
matrix in Keldysh space. Equation (^) generalizes the Landau-Boltzmann transport equation for normal Fermi liquids 
to the superconducting state. In our model we have a transport equation for each layer, specified by the discrete 
index £. The spectrum and distribution of quasiparticles in layer I are affected by the in-plane vector potential, 
Ai(t), the scalar potential, &e{t), the pairing field, Ai(s,t), and the scattering self-energy, te(s;e,t). The coupling of 
quasiparticles (of charge e) to an electromagnetic field appears explicitly in equation (|3j), and also implicitly through 
the field dependence of the self-energy, S£. The Fermi velocity, Vf, is a two-dimensional vector parallel to the planes. 
The commutator in equation (||) is defined as [a, b)® = a ® b — b <g) a, where the symbol <g> denotes a folding product in 
the energy-time domain (see Appendix), together with matrix multiplication of Nambu-Keldysh matrices. For details 
of this formalism see Refs. ^3|-|37|. 

In the weak-coupling approximation, the retarded, advanced, and Keldysh components of the pairing field take the 



form 



Af A ( S ;i)= J ^-.(v(s, S >)f e K (s';e,t)) si , (5) 
Af(s;t) = 0. (6) 
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Here V(s, s') is the pairing interaction, which determines both T c and the symmetry of the order parameter, and f K 
denotes the off-diagonal part of the Keldysh propagator, g K . The notation (F(s)) s denotes a Fermi surface average, 
/ cPsn(s)F(s), where n(s) describes the local density of quasiparticle states on the Fermi surface, normalized so that 
J d 2 s n(s) = 1. For the case of isotropic s-wave pairing, the pairing interaction takes the form V (s, s') — Vq, while for 
a d-wave model, V (s, s') favors a pairing field with d x 2_ y 2 symmetry, although essentially all of the results presented 
here are valid for any unconventional order parameter with line nodes on the two-dimensional Fermi surface. 

The scattering self-energy of layer £, Eg = s| + + , describes both in-plane scattering ) and interplane 
scattering (p e t±x)- In-plane scattering processes are taken into account via a self-energy of the form 

flj (s; e, t) = Ci it (s, s; e, t) , (7) 

where Cj is an effective concentration of scattering centers. We consider isotropic in-plane scattering, in which case 
Uo = uq 1 is an isotropic scattering potential. The s-dependence drops out of the in-plane t matrix so that 

k (e, t) = uo + u <g> [N f (g t (s; e, t))J ® i e (e, t) , (8) 

where Nf is the density of states at the Fermi energy. The scattering self-energy, sf, then contains only two inde- 
pendent parameters. Following Buchholtz and Zwicknagl,E£l we choose these two parameters to be the normal-state 
scattering rate, 

h_ _ 2n Cl N f ul ^ 



Tl 



and the normalized scattering cross section, 



]| 1 + (NjuonY 



{NfUo7T)2 2 , (10) 



1 + (NfUoirY 

which is the true cross section divided by the cross section in the unitarity limit, a^J = 4/fcy. Thus, the normalized 
cross section ranges from a — for weak scattering (Born limit), to a = 1 for strong scattering (unitarity limit). We 
have formulated the in-plane scattering self-energy in terms of a t matrix. In this formulation, the self-energy covers 
impurity scattering, including resonant scattering, as well as any other type of scattering process, elastic or inelastic, 
as long as it can be approximated by an effective lifetime and an effective cross section. For instance, electron-phonon 
scattering in this approximation corresponds to er = (Migdal's theorem) and a temperature dependent lifetime, 
r,|(T). 

Since we neglect coherent transport along the c-axis (i.e., we set the Fermi velocity along the c-axis to zero), 

the interlayer scattering self-energies, T, e are the only source of interlayer coupling in the model. We assume this 
coupling to be weak (nonresonant), and thus describe interlayer scattering by a self-energy in the Born approximation, 

£f t t ±1 (s;e,t) = U iii±1 (t)® A^—L—^^je^ ®^ ±1 (i). (11) 

The interlayer self-energy is parametrized bv-an effective lifetime, t_l, which is the characteristic time between con- 
secutive scattering processes between layers.eSl The gauge operators, Ui/±i, and interlayer vector potentials, A z ll±l , 
are defined by 

Ue,e±i(t) = exp ( -<|^A|, ±1 (i)fa ) , (12) 
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he 
±i) 



^W*) = 3/ dzA*(z,t), (13) 

« J id 

where d is the layer spacing. The effective interlayer scattering lifetime, t±(s. s'), is generally anisotropic. We describe 
this anisotropy by 

I _ i exp [7 cos(p(s) - (pis'))] - 14 . 



t±{s,s') t± ((exp [7 cos (<p(s) - (pis'))]) s ) s > 

The real space angles ip(s) and (pis') give the in-plane directions of the quasiparticle Fermi velocity before and after 
scattering to an adjacent layer. The parameter 7 specifies to what degree the scattered electrons "remember" their 
initial momentum. Isotropic scattering corresponds to 7 = 0, while extreme forward scattering corresponds to 7 — > 00. 
The scattering rate is illustrated in figure ([[]) for several values of the parameter 7. Note that in the interlayer diffusion 
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model one needs.-an anisotropic c-axis scattering rate to get a nonzero interlayer Josephson coupling_u/ith a d-wave 
order parameter O However, higher order tunneling processes also give a nonzero Josephson couplingo 

Here we are concerned with the current response to a weak electric field for a layered d-wave superconductor. The 



in-plane current is carr ( 
of quasiclassical theor; 



quasiparticles moving with an in-plane velocity V/, and is given by standard equations 
de 



u(t) = eN f I —( Vf ( s )Tr[f 3 gf(s;e,t)]) s 



(15) 



The interlayer current, on the other hand, is of a different nature. It originates from random scattering processes of 
quasiparticles between adjacent layers. The microscopic expression for the interlayer current density was derived for 
isotropic interlayer scattering in Ref. |^, and is generalized below to anisotropic scattering, 




K \ 



e+i ® 9i 



e+i ® 9i 



9 1 ® ^t,i+\ 



gf )}} [ 



(16) 



Equations © and © cover the current response for both normal and superconducting states, in equilibrium and 
nonequilibrium situations. 



B. Equilibrium Solution 



In the absence of external perturbations, the quasiclassical transport equation and normalization condition for the 
equilibrium retarded and advanced propagators become 

•R.A 



eh - A R > A (s) 



(s;e),g RA ( S ;e) 



= 0, 



(17) 



[g RA {s^] 2 = -^i, 



(18) 



where the ® products reduce to simple matrix multiplication. We dropped the layer index £, since all layers are 
equivalent in equilibrium. In the real gauge, the pairing field has the simple form, A R ' A (s) — iAo(s)r2. Choosing this 
gauge, we solve equations (0) and Jig), and obtain 



(s;e) = -7T- 



~R,A 



( S ;e)f 3 -iA*- A (a;e)f 2 
A R ' A {s;e) 2 -e R ' A (s;e) 2 



(19) 



where the e and A can be interpreted as the "renormalized" excitation energy and gap function due to scattering 

~ R A 

processes. The renormalized quantities are determined by the scattering self-energy, £ : , and are given by 



~e R > A (s;e) =e-^Tr 



1 



* -R,A , ,' 
r 3 S (s;e) 



A R > A (s;e) = A (s)--Tv 



~ ~R,A , \ 

r 2 s (s;e) 



The in-plane scattering self-energy can be evaluated by solving equation (|^) to obtain 

(g R > A (s;e)) s 



~R,A , \ 



2th 



1 — o 



l + (ir-Hg R > A (s;e)) s y 



- 1 



while the interplane self-energy has the simple form 



-R.A, , 



7r\rj_(s,s') /, 



(20) 
(21) 

(22) 
(23) 



From the solutions for the equilibrium retarded and advanced propagators and self-energies, one can calculate the 
equilibrium Keldysh components, 



g (s; e) = tanh 



T, K (s; e) = tanh 



2k B T 
e 

2fc R T 



) (g R (s;e)-g A (s;e)), 
) (t R (s;e)-t A (s;e)) 



(24) 
(25) 
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C. The d-wave Model 



So far the presentation of the quasiclassical theory has been fairly general. To make further progress we must 
specify the shape of the Fermi surface and the form of the pairing interaction, V(s, s'). We consider a simple model 
that incorporates d-wave pairing. We assume a cylindrical Fermi surface (in which case the Fermi surface position, s, 
is replaced by the azimuthal angle, 0), and a pairing interaction of the form 



V (0, <p') = V cos (20) cos (20'). 



(26) 



Note that the basis functions for this pairing interaction possess the usual fc 2 — fc 2 anisotropy. The angular dependences 
of the self-energies are now determined, and the remaining task is to solve for the magnitudes and energy dependences. 
For a system in equilibrium, and in the absence of any perturbations, we find from equation (0) that the pairing field 
has the form 



iA cos (20)f 2 , 



and from equation (22) the in- plane scattering self-energy becomes 



*r,a , % 
2 II ( £ ) 



1 



2ri, 



where 



(gi' A (fce))<t> = 



-Tre R ' A 



(e) 



2- 



A 2 cos 2 (20) - e R ' A (e) 2 



(27) 



(28) 



(29) 



Note that £> R ' A (e) has no t 2 component because the Fermi surface average of A"' A (0) is zero. However the r 2 
term does not necessarily drop out of the interplane scattering self-energy because of the anisotropy of the interplane 
scattering rate, t±(s,s')~ 1 . 

In the calculations that follow we assume that the interplane scattering rate, 1/tj_, is small compared to the in- 
plane scattering rate, in which case we neglect corrections of the order 1/t_l in the renormalized self-energies and 
retain only the leading contribution in 1/t± that determines the c-axis current response. The equilibrium propagator 
can then be written as 



9 



R,A , 



0; e) 



~R,A 



(e) f 3 - iA cos (2^) f 2 



A 2 cos 2 (20) 



?R,A 



(30) 



where the renormalized excitation energy is given in equation (pOj). Normally one eliminates the pairing interaction, Vq, 
and frequency cutoff in favor of T c , which is then a material parameter taken from experiment. The gap parameter, 

- R A 

Ao(T), is then calculated self-consistently in terms of T c and the parameters defining the self-energies Em ' (e). 
In the absence of detailed information on the dominant inelastic processes we have effectively made relaxation time 
approximations for (s; e), which are parametrized by temperature-dependent relaxation rates, r^ 1 and rj^s, s'). 
Thus, we take rn (T), t±(T), and Aq(T) as material parameters and calculate measurable properties below T c in terms 
of these parameters and the assumed pairing symmetry. Given a gap Aq, one must solve equations (|2^) and ( |30| ) 
self-consistently for £ R ' A (e). The equations can be solved numerically via an iterative technique to obtain the density 
of states 



AT, 

N(e) = ^(gi^e))^ 



(31) 



The equilibrium solutions for i R ' A (e) and g R ' A (0; e) are also needed to carry out the current response calculations 
described below. Although the density of states, for an unconventional superconductor with lines of nodes, has been 
known for some time (see review [l2|) we find it useful to discuss our results for the conductivity in conjunction with 
the density of states. 



III. LINEAR CURRENT RESPONSE 



Calculation of the infrared conductivity requires an appropriate formulation of linear response theory. The quasi- 
classical formulation of linear response has been developed by several authors (cf. Ref . [43]) . We give a brief summary 
of an efficient formulation of linear response theory that utilizes tke normalization conditions and related algebraic 
identities to solve the linearized quasiclassical transport equations S3 We consider a time-dependent electric field that 
is uniform in space. This is a good approximation for high-T c superconductors with a penetration depth much larger 
than the mean free path and the coherence length. 
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A. In-plane Conductivity 



For a uniform field it is convenient to write equation (||) in the compact form 

[eh-SH-t,g] 9 = 0, (32) 

where Sv is the perturbation, and the pairing field, A, has been included in the full self-energy, S. For the electrical 
conductivity the perturbation is of the form 

5v R < A ( S) t) = -~v f (s)-A(t)T 3 , (33) 
c 

Sv K (s,t) = 0. (34) 

We assume that the equilibrium solution, denoted by So and go, is known, then linearize the quasiclassical transport 
equations in terms of the perturbation, 8v(s;t), and the response, Sg(s;e,t). In general these linearized transport 
equations involve a self-energy response, St, (s; e,t), corresponding to the vertex corrections in a Green function 
formalism. However, for the case of uniform perturbations and isotropic scattering rates, the vertex corrections to the 
transverse conductivity vanish. Note that in calculating the in-plane current response, we neglect the contribution of 
the interlayer scattering self-energy to So. We write the matrix propagator in the form 

9 (s; e, t) = g (s; e) + Sg (s; e, t) , (35) 

then expand equation ( |32"|) to first order to obtain the linearized transport equation, 

[ef 3 - S (s; e)] ®8g (s; e, t) - Sg (s; e, i)<g> [ef 3 - S (s; e)] 
= 5v (s; t)®g a (s; e) - g (s; e)®6v (s; t) . 

(36) 

In addition, we expand the normalization condition through first order and obtain 

go(s;e)®go(s;e) = -7T 2 i, (37) 
50 (s; e) ® Sg (s; e, t) + Sg (s; e, t) ® g (s; e) = 0. (38) 

The (g> product denotes a folding product in the energy-time domain, together with matrix multiplication of Nambu- 
Keldysh matrices. Fourier transforming the time variable turns the energy-time folding in equations ( (36|) and fl38|) 
into a simple product, with an energy shift of magnitude 4p in the equilibrium quantities. This simplification is a 
consequence of the time independence of the unperturbed propagator and self-energy. Thus, equations ( |36| ) and (|3^) 
become 

[(e+^)f 3 -S (s;e+ff)] ®Sg(s;e,uj)- Sg(s;e,u>)® [(e- f 3 - S (s;e- %i)] 

= <Jfi (s; w) ® g {s; e - ^) - .90 («; c + 4f ) ® <5v (s; w) , (39) 



50 (s; ®Sg (s; e, w) + Sg (s; e, w)® 5o (*; e- fl) = 



(40) 



where the ® product in equations (p7|), (J39|), and ( pOj ) should now be interpreted as a simple matrix-product of 
Nambu-Keldysh matrices. 

Equation (J39]), supplemented by equations d37j ) and (^), amounts to 16 coupled linear equations for the 12 inde- 
pendent components of the Nambu-Keldysh matrix Sg. These components are conveniently described by the thrpft 
2x2-Nambu matrices Sg R , 5g A , and 5g K . The linear equations can be solved either directly, via matrix inversion,c2l 
or, alternatively, by using special algebraic identities connecting the propagators and self-energies, together with the 
normalization conditions (37) and (BQ)E3 The identities are based on the following relation between the terms in the 
square brackets in equation (p9|), and the unperturbed propagators, 



£T3 



C 



R,A 



(«; e ) % ' A ( s ; e ) + D 



R,A 



(s;e)l, 



(41) 



where C R ' A (s; e) and D R ' A (s; e) are c-number functions that are uniquely determined by the solutions for So and go- 
The general solution of the linearized transport equation can be written in terms of the unperturbed propagators, the 
perturbation, and the functions C R ' A , and D R ' A , 



~R,A 



(s;e,u) = - 



C+ {s;e,uj)g ' (s;e+— ) 



t^R.A, \ 

Dj (s;e,w) 



7r^G + ' (s;e,uj) + Dj (s;e,w) 



Sv (s; w) % 



-R,A 



hu\ ~R,A 1 



■^)Sv(s 



(42) 
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and 



Sg K (s; e, w) = tanh (^Jr J (*5 e, w) - tanh \^Jrj Sg A (a; e, w) 

+ ( tanh (^)" tanh (^))^ (S;£ ' W) ' ^ 

where we have written the. Keldysh propagator, 5g K , in terms of Sg R ' A , and the "anomalous" response function 
introduced by Eliashbergj23 which has the solution, 

sg a («; e, «) = - ^'^f^T^^T^ ^ (s . u) § a {s . e _ ¥) _ g n (s; e+¥) 5€ (s . w)] (44) 

7r 2 C7£ (s; e, u/) + £> £ L(s;e,w) " 

The functions C+ ,j4 ' a , D R ' A ' a in equations ( f42| ) and are defined in terms of C fl,A and D^'^ 1 of equation ( |4l| ) as 
follows, 

C^(*; e , W ) = C^( S ;e+4f)+C^( S ;e-4f), (45) 
D*' A (s;e,u/) = D fl ' A (s;e+^)- J D i *» A ( S ;e-Af), (46) 
C£(*;e,w) =C R (s;e+!f)+C A ( S ;e-!f), (47) 
££(«;e,w) = £ fl (s;e+^f) -£ A (s;e-^f) . (48) 

Of central importance is the Keldysh propagator, Sg K (s; e,u>), since it gives directly the linear response of ob- 
servable quantities, such as the charge density and the current density, to the perturbation, Sv. Inserting Sg K from 
equation (|4^), together with Sv from equation (|33|), into equation ( |l5| ) for the current gives the in-plane conductivity 

2N f e 2 v 2 f f de f d<f> 2/1 .\ , C£(0;e,w) 
all w = — L L — -f. cos 2 U)< tanh - 2 + , ' ? 



tanh 



x (0;e-4f ) S * (0; e+Af) + /<f (0; e-^) /<f (0 



e+%0+7r 2 l 



2fc B T y ^C*(ct>;e,ujy +D A (cf>;e,u>y 

x [ 5o A (0; e- If) g A (0; £ + 4f ) + f A (0; e- 4f ) f A (0; e+ ^ 
'e+^\ _ tanh / e_^\ \ C%{^) 
v 2fc B T j ^ 2fc B T ^Ci (0; e, a;) 2 + D°_ U; e, W ) 2 



x [<tf(0; C -^) S *(0; e +M)+# (0; e -^) /« (0; e + If ) + tt 2 ] . (49) 



where ffjf'^ and denote the upper-left and upper-right elements of the equilibrium Nambu matrices g^' , re- 



spectively. Results for the in-plane conductivity of a <i-wave superconductor are discussed in section [V 



B. Interplane Conductivity 



Now consider the current response for a weak electric field polarized along the c-axis, E z — — -d t A z {f). From 
equation {UM, the c-axis current response becomes 



5f(t) = — 7 



eiVfd /" de 



t Tr f 3 



(5s (s;e,i) ,go(s;e) 



(50) 



where we neglect the response of g to A z . The interlayer self-energy is evaluated by expanding the gauge matrices 
defined in equation (12) to first order in the perturbation, 



ied I 1 \ 

5± ± (s; e ,i) = - — ( — [A z (t)r 3 ® § (s'; e) - g (s'; e) ® A z (t)f 3 ] ) 

2ttc \tx(s,s'J / s , 



(51) 
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Substituting equation (|T]) into equation (50), and Fourier transforming the time variable, we obtain the linear response 
of the c-axis current, 5j z {uj). The expression for 5j z (u) can be formally separated into a quasiparticle term, Sj z p (u>), 
and a Cooper-pair term, 6j*(w), or equivalently, two contributions to the interplane conductivity, a = o^ p + er^, 

_X /. .\ N f eW f de [dcf> [d<f>> 1 f r R n A h 1 

x [ 9 K(4,>,e+!f)-gK(4>',e-!f)] 
+^(0, £ +%i) [^(0', e +^)-^(0' !e -^)] -sf(0, e -^) [ ff( «(0', e + ^)- 5o V, e -^)] }, (52) 

and 

, Nff£_ [ de f d4> f dfi 1 f ■■ fl ft A ^ , 

x[/ ^ e + 4f) + / ( f(<//,e-4f)] 
+/«f (0,6+%i) [/oV^+^O + ZoV,*-^)] +/<f(<^-^) [/^'.e+^O+ZW^-^)] }• (53) 

Note that for an isotropic scattering rate (rj 1 (0, 0') = tJ 1 ) the quasiparticle contribution to a 1 - is finite, while the 
Cooper-pair contribution to the conductivity vanishes in the case of a d-wave gap function. This happens because the 
angular average of the d-wave gap-function, and therefore (f^' A ' K )<f,, vanishes by symmetry. Since the zero-frequency 
limit of a^ p determines the dc Josephson current, anisotropic interplane scattering is necessary to obtain a nonzero 
Josephson coupling.cHl 

The low- frequency limit of the imaginary part of er (T, ui) determines the London penetration depth, Aj_(T), 

' 4 " limujImo~ L (T,w). (54) 



Ai(T) c 2 

For s-wave pairing, we obtain a relation between the normal-state c-axis conductivity, energy gap and c-axis penetra- 
tion depth which is equivalent to that derived for superconducting alloysJU 

1 =^(T)Ao(T)tanhf^V (55) 



Aj_(T) he 2 ny J >y ' V 2fc B T 

For incoherently coupled layers cr^{T) =2e 2 Nfd 2 /rj_(T). This relation is useful for estimating the London penetration 
depth from measurements of the normal-state dc resistivity and T c , and it provides an important consistency check 
on the interlayer tunneling model with s wave pairing. 



IV. RESULTS 

The key complication in almost any quasiclassical calculation is the requirement that g and £ must be calculated 
sclf-consistently. While this is a formidable impediment to analytic calculations, one can easily solve the self-consistent 
equations using elementary numerical techniques with relatively limited computing power. The results presented here 
were obtained from numerical calculations carried out on desktop workstations. 



A. Model Parameters 



Our model for the cuprate superconductors contains, as a minimal set, five phenomenological material parameters: 
the transition temperature, T c , the total density of states at the Fermi level, Nf, the Fermi velocity, Vf, and the 
in-plane and interplane scattering lifetimes, rii and t± for the normal state. All of these quantities can be deduced 
from normal-state measurements. In order to accommodate d-wave pairing in the interlayer scattering model, we 
introduce two additional material parameters: the normalized scattering cross section, a, and the interlayer scattering 
anisotropy parameter, 7 (see equation (|l4j). Our goal is to check the key assumptions of the interlayer diffusion 
model; (i) that the superconducting planes are Fermi liquids, and (ii) that the interlayer transport is dominated by 
incoherent scattering processes. We use normal state data as input to calculations of superconducting properties, 
which then allows us to check the consistency of our model. We examine models based on both s-wave and d-wave 
order parameters. 
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The best characterized high-T c material to date is Y-Ba-Cu-O. The procedure for obtaining approximate values for 
the material parameters (excluding a and 7) was discussed in detail in Ref. [| Using the superconducting transition 
temperature, T c — 92K, the Sommerfeld constant, 7 S = 200/j J/{K 2 cm 3 ), the in-plane resistivity, Q\\(T C ) = 50/if2cm, 
the normal-state Drude plasma frequency, hui p — 1.5eV, and the resistivity ratio, [Q±/q\\) t = 200, we obtain the 
material parameters given in table |[ We can then calculate superconducting properties within our model, then check 
for consistency with published experimental values. Here we extend the treatment of Rcf. ^ to include d-wave pairing. 

Our results, for equilibrium supercurrents, are given in table |n|. The lower bounds are estimates obtained assuming 
a vanishing in-plane scattering rate, while the upper bounds were obtained using the lifetimes given in tabic |[ and 
a = 1 for the d-wave model. The c-axis quantities were calculated using a moderate anisotropy parameter of 7 — 0.5. 
The "in-plane" quantities (A|| and dH^/dT) agree well with experimental values for either an s-wave or a d-wave gap, 

while quantities involving c-axis transport (Aj_ and dH^/dT) yield better agreement for s-wave pairing. We interpret 
the in-plane agreement as support of the hypothesis that the superconducting layers of the cuprate compounds are 
well described by a Fermi liquid theory (regardless of the symmetry of the order parameter). The poor agreement 
obtained for d-wave pairing, when c-axis properties were considered, can be interpreted as a failure of either the model 
of incoherently coupled layers, or the model of pure d- wave pairing, or both. 

In addition to Y-Ba-Cu-O, we have compare d . the i nterlayer diffusion model with recent data on Ar-annealed Bi-Sr- 
Ca-Cu-0 and 02-annealed Pb-Bi-Sr-Ca-Cu-0;E3E3e3 penetration depth, Aj_, critical voltage, V c , and critical current, 
I c , measured at A.2K . The data is insufficient to carry out a complete analysis, but we can check the consistency of 
the theory with regard to the c-axis properties. Two independent relations can be derived for the c-axis transport 
properties at low temperatures, T <C T c , 

^Mj^d=l, (56) 



Air 2 A^Aq _ f 1, for s-wave pairing, 

(57) 



he 2 q£ \ Rd, for d-wave pairing, 

where d is the layer spacing, and Rd is a dimensionless number. 

Equation ( [56] ) relates the interlayer Josephson critical current to the penetration depth and is independent of the 
pairing symmetry. The second relation ( |57|) is sensitive to the symmetry of the order parameter and depends critically 
on the interplane anisotropy parameter, 7, in the case of d-wave pairing. In the limit of isotropic scattering (7 — > 0) 
Rd — > 00; however, Rd is hnitc for anisotropic scattering (7 ^ 0), and can be ~ 1 fan extreme forward scattering 
(7 — ► 00). Combining equations (|5^) and (^) yields the Ambegaokar-Baratoff relatione!} 

2e 1 j, / 1, for s-wave pairing, 



7rAo \ Rd ' for d-wave pairing, 

which is independent of 7 for s-wave pairing, but not for the d-wave model. Taking values for Aj_, j^r, and w 
V^~/j^d, from the experimental data, and estimating Ao from T c and the clean- limit weak-coupling relations; — 
1.76k B T c , A { Q d) = 2.Uk B T c , we can evaluate the ratios given in equations (E^)-(E^) for the Bi-Sr-Ca-Cu-0 compounds. 



The results are summarized in table III. The data on Ar-annealed Bi-Sr-Ca-Cu-0 are consistent with s-wave pairing, 
and also with d-wave pairing provided the interplane tunneling is extremely anisotropic (7 > 5). Neither the s-wave 
nor d-wave model accounts for the 02-annealed data; thus, the interlayer transport mechanism in these compounds is 
not well described by incoherent tunneling—. Such systems may be better described bypXesiipaxLt tunneling via impurity 
states localized between the Cu02 layers,cZl or perhaps other transfer mechanisms. Ililj'EjEfluJ 

The analysis of the equilibrium properties implies that Fermi liquid theory provides a sound basis for the description 
of the in-plane properties of high-T c superconductors, independent of the symmetry of the superconducting order 
parameter. The interplane properties, on the other hand, depend strongly on the gap symmetry. The experimental 
situation is also not simple; results differ depending on the compound and material preparation. For some samples we 
find reasonable agreement with the universal relations implied by the interlayer diffusion model, while other samples 
display significant deviations. These discrepancies are present for both s-wave and d-wave models, and suggest that 
the model of incoherently coupled layers is an oversimplification of the c-axis charge transfer mechanism in certain 
high-T c materials. Alternatives are, for instance, resonant coupling, coupling by pair tunneling, coherent coupling, or 
a combination of different types of coupling mechanisms. 



B. Conductivity 



The procedure used to solve equations ( P0|) and (|3y) for e R,A , which then yields the equilibrium propagator and 
self-energy, was described at the end of section II. In solving for the equilibrium solutions we neglect the interlayer 
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self-energy term. This is justified because the interlayer scattering rate is much smaller than the in-plane scattering 
rate. To calculate the in-plane conductivity, the remaining step is to carry out the integrations in equation ([l9|). A 
key feature of the d-w&ve model is the presence of low-energy excitations near the nodes of the gap. The phase space 
for the d-wetve model is similar to that of the three-dimensional polar (p-wave) model with a line of nodes. We note 
that there is a good qualitative agreement between the infrared conductivity calculated for a three-dimensional polar 
state.Ej'EZI and the two-dimensional d-wave state considered here. Both states possess lines of nodes in the gap. 

In figure (^a) we show the density of states for a d-wave order parameter in a fairly clean sample (i.e., a = 
h/2irrAo = 0.01), for several values of the .normalized cross section. The key feature to note is the enhancement 
of the density of states at low frequencies. He3c3 This feature is strongest in the unitarity limit, and softens as one 
goes to lower values of <t, until it vanishes completely in the Born limit. This enhancement is interpreted as being 
due to bound-states, localized at the scattering centers, whose origin lies in local pair-breaking effects.EIH At a finite 
concentration of scattering centers, the bound states overlap and form "bound bands" . 

The bound states are optically active; the matrix elements coupling the impurity bound states do not vanish as 
q — > 0. This is evident from figure (|^b), where we show the in-plane conductivity for a sample with the same 
parameters as in figure |^a), at T = 0. In the unitarity limit [a — > 1) the absorption is substantially enhanced at 
low frequencies compared to the Born limit (a — > 0). The enhancement is due to transitions within the bound band. 
Note that the enhancement of the conductivity at low frequencies, which shows the strong optical activity of the 
bound states, is much more pronounced than the enhancement in the density of states. Furthermore, the in-plane 
conductivity approaches the universal limit obtained by LeeJE3 lim^^o a 

II (w) ~ N f e 2 v 2 f h/7rA , independent of the 
scattering cross section (er), as shown in figure (^|b), and scattering rate (a). In addition to the strong absorption at 
low frequencies there is also a small feature at Tiuj ss Ao , which corresponds to the excitation of a quasiparticle from 
a bound state to the continuum states at energies Tiuj ~ Ao- Finally, we note that there is no discernible feature at 
Tiuj = 2Ao, in contrast to the case of s-wave pairing. Thus, optical experiments are not expected to provide a good 
means of measuring the maximum gap in d-wave superconductors. 

In figures ([|a) and (|3|b) we show the density of states and in-plane conductivity for a sample with a larger scattering 
rate (i.e., a — Tl/2tttAq — 0.1). In this case the bound states have been merged with the continuum states and are no 
longer distinctly visible in the density of states. However the effects of the bound states can still be seen in the infrared 
conductivity. Note the different scale between fig. (J^b) and fig. (||b) is due to the change in the normal-state Drude 
conductivity. The universal value of the conductivity for the superconducting state implies lim^^o c" (w)/crjj(0) ~ 2a. 
Note also that the absorption at low frequencies is still much larger for resonant scattering (tr — ► 1) than it is in the 
Born limit (a = 0). 

In figure (Ji|a) we show the temperature dependence of the in-plane conductivity for the unitarity limit and a = 
Tl/2tttAo = 0.01 as in figure (0). Typically the scattering rate, r _1 , will itself depend on temperature. However, this 
temperature dependence, which is obtained from experimental measurements in our model, has been neglected here. 
Note that the gap-like feature at Tiuj « Ao is strongest at low temperatures, and can be as much as 50% higher than 
the normal-state Drude value of the conductivity. At higher temperatures we find an upturn in the conductivity at 
low frequencies, which yields a coherence peak (i.e., the conductivity exceeds the normal state value). 

The c-axis transport calculations were carried out with a moderate value of the scattering anisotropy parameter, 
7 = 0.5 (see equation (|l4|), unless otherwise stated. In figures (^c) and (||c) we show results for the real part of the 
c-axis conductivity at T = 0, for two values of the scattering lifetime, and for both s-wave and d-wave pairing. Note 
the onset of dissipation for isotropic s-wave pairing at Tiuj = 2Aq. It is interesting to note that although the dynamical 
tunneling currents, for the case of an s-wave gap, obey the frequency dependence of "dirty" superconductors, their 
dependence on the scattering lifetime is just the inverse (i.e., a 1 - oc tJ^ 1 ). Thus, in the interlayer diffusion model, 
an increase in interlayer scattering increases the c-axis charge transfer (regardless of the order parameter symmetry) , 
whereas it reduces the charge transfer in "dirty" alloys. It is also significant that the c-axis conductivity in the 
interlayer diffusion model does not tend to a universal value in the limit uj — > 0. 

A d-wave superconductor has residual dissipation even at u> — 0, due to low energy excitations near the nodes in the 
gap, which is significantly enhanced in the limit of strong scattering. The knee that develops in the limit of resonant 
in-plane scattering (er — > 1) at Tiuj Ao is due to empty bound states, as discussed for the in-plane conductivity. 
Thus, it may be possible to shed light on the strength of in-plane scattering by detecting this common feature in both 
conductivities. In figure (Qb) we show the temperature behavior of the c-axis conductivity for the unitarity limit, 
and a — Ti/2ittAq = 0.01. It saturates at zero-frequency and has a washed-out minimum, which shifts slightly with 
increasing temperature to higher frequencies. 

C. Reactance and Low- frequency Plasma mode 

The imaginary part of the conductivity (reactance) for an s-wave gap is well described by a term of the form 
cr^-(w) ~ A/ttuj for < Tiuj < 2Ao, where A is determined by the dc penetration depth, \f~A — c/2A^. This 
quasihydrodynamic approximation for the conductivity is valid at sufficiently low frequency, when quasiparticle pair 
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production is negligible and the current response is dominated by the London supercurrent. For an s-wave gap at 
T — pair production is absent for htv < 2Ao, so the London limit is a good approximation even for hu> ~ Ao. 
However, the reactance for a d-wave gap deviates from the w _1 behavior at frequencies much smaller than the gap, 
dropping much faster than the s-wave result, and crossing zero at Tilo ~ 0.3Ao (see figure (H). The breakdown 
of the London limit develops at much lower frequencies for the d-wave gap because quasiparticle pair production 
is present at any frequency for a gap with nodes. The fraction of the Fermi surface for which pair production is 
energetically allowed is 4> u ~ {hoj /tv Aq) . Pair production leads to a reduction in the supercurrent, i.e., a softening 
of the condensate, and of course to dissipation at all frequencies below the maximum gap of 2Ao. The softening 
of the condensate current response appears as a reduction in the Cooper pair contribution, to the conductivity 
compared to the quasiparticle part, a^ p . The dashed and dotted lines in figure (^) correspond to the condensate and 
quasiparticle contributions to the reactance, respectively. For the case of an s-wave order parameter the condensate 
contribution has the highest weight at least through frequencies on the order of the gap amplitude, while the reactance 
for a d-wave superconductor is dominated by quasiparticle excitations at frequencies much smaller than the gap. 

This softening of the condensate response at frequencies below ~ Ao, by the quasiparticle excitations, has important 
implications for the low- frequency plasma mode in strongly anisotropic superconductors. The large anisotropy of 
the penetration depth in the cuprate superconductors (Aj_/A| 3> 1) leads to a low-lying plasma mode along the 

c-axis with frequency uj pc < lo p A||/Aj_,Ej where io v is the in-plane Drude plasma frequency. The anisotropy of 
the cuprates is sufficiently large so that the collective mode lies below 2Ao.E3 The transverse dielectric constant, 
e±(u>) = £oo + Ania^ {id) / lo , vanishes at u> = uj pc ; thus, for low damping the plasma frequency is obtained from the 
root of the real part of s±(llj). The frequency, u>*, at which Imo" L (w) crosses zero provides an upper bound for the 
plasma frequency, u> pc < lo* . Shown in the insert of fig. (^|) is a plot of Im a J ~ (u>) for both the s-wave and d-wave models 
with the same zero-temperature penetration depth, Aj_(0) (that is, we scaled Imcr^(w) of the d-wave model so that 

Imcn^(O) = Imcr^(O)). Since the root, u>*, of the reactance is maximal for the clean limit we only display the result 
for a — * 0. In the d-wave case we can deduce an upper bound for uj pc from thpJutjcrccpt of the reactance with the 
w-axis, which for 7 = 0.5 falls below the range of experimentally observed valuesaii] (denoted by the shaded region). 
Note, that uj pc increases with increasing 7, but is always below the s-wave result. We conclude from our analysis of 
the reactance that the London approximation for the c-axis conductivity (and dielectric function e(w)) breaks down 
in the case of an incoherently coupled stack of two-dimensional d-wave superconductors, if the interlayer scattering 
is moderately anisotropic (7 < 0.5). Finally, we note that the pair production and softening of the condensate at 
relatively low frequencies leads to an anpmalous temperature behavior of the surface resistance, which increases with 
decreasing temperatures well below Ao.E3 

V. CONCLUSION 

We have calculated the in-plane and c-axis infrared conductivities for a layered d-wave superconductor within the 
Fermi liquid theory of superconductivity, including the effects of in-plane and interplane scattering. An important 
effect on in-plane scattering, for d-wave superconductors, is the formation of quasiparticle bound states that are 
localized at the scattering centers. These bound states are clearly visible in the density of states, but they have 
a far greater effect on the infrared conductivity (both in-plane and c-axis) because the localization of these states 
at impurities leads to a large optical activity. The effects of the bound states are most pronounced, for resonant 
scattering centers; Zn impurities in Y-Ba-Cu-0 are argued to be candidates for resonant scatterers.cMj Observation 
of the conductivity from impurity states in the superconducting state of high-T c materials could provide important 
evidence for d-wave pairing. 

Our calculations of the interplane conductivity also show that the transport properties of a layered d-wave super- 
conductor are very sensitive to the nature of the interlayer coupling. For instance, incoherent interlayer coupling, 
mediated through isotropic scattering processes, yields a vanishing Josephson coupling because the angular average of 
the d-wave gap function vanishes for symmetry reasons; however, for highly anisotropic scattering (conserved in-plane 
momentum) the Josephson coupling becomes of the order of the s-wave value. As a result, it is difficult to reconcile 
the model of incoherently coupled layers, together with d-wave pairing, with the currently available c-axis transport 
data. However, s-wave pairing in the interlayer diffusion model often yields a reasonable fit to the c-axis data. Thus, 
a better understanding of the microscopic nature of the interplane transport mechanism is needed in order to use 
c-axis transport measurements as probes for the order parameter symmetry in high-T c superconductors. 
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APPENDIX: 



The study of nonequilibrium properties of a Fermi liquid, requires knowledge of the dynamics of both the distribution 
of low-energy excitations, as well as the spectrum of these states. This may be accomplished within a framework 
introduced by KeldyshBa This formulation allows us to treat the spectrum and distribution function on the same 
footing, and to obtain the dynamics from the standard array of Green's function methods. We provide a brief summary 
of the definitions for the Green's functions used in this paper, as well as the microscopic definitions of the quasiclassical 
propagators. 

The two-component Nambu-field operator in layer I is defined as ^] (x, t) = f^lf ( x > > V^j. ( x : f° r s P m singlet 
superconductors. The retarded, advanced, and Keldysh Green's functions are defined as follows ill 

Gf (x, t; x', t') = - ~0{t - t') ( { 9 t (x, t) , 9\ (x ', t>) } ) , 
G A (x,t;x',t') = i*(f - t) (x,i) (x ',*')}) , 

Gf (x,i;x',0 =-i([* / (x,t),*J(x , ,t')])- (Al) 

The "hat" denotes a 2 x 2 matrix in particle-hole (Nambu) space. These propagators contain information on the phys- 
ical quantities of interest. The diagonal component of Gf' A ' K describes the spectrum and distribution of quasiparticle 
excitations, while its off-diagonal component, F^' A ' K , describes the pairing correlations. 

The quasiclassical propagators, gf' A ' K , are obtained from the microscopic propagators, Gf'' K , by integrating-out 
the short wavelength, high-energy information. In the mixed (p, R) representation, the formal definition of g^' A ' K is 

~R,A,K / t> ,\ 1 / ,> « AM,A,K 



(s, R;e,t) = - d£ p r 3 G?' A ' n (p, R; e, t) , (A2) 

a J -lu c 

where a is the quasiparticle renormalization factor .H In this paper we are interested in spatially homogeneous systems 
in which the I and R dependence drops out of equation ( |A2| ). Thus, our propagators are functions of only the Fermi 
surface position s, the quasiparticle excitation energy, e, and time t. In the transformation to the mixed representation, 
the Green's function products appearing in the Dyson equation become energy-time folding products which we denote 
by the symbol ®. A useful formal representation of this folding product is: 

a(s;e,f)<8 b(s;e,t) = e* (0 ° a * - d ? d ^a(s; e, t)b(s; e, t), (A3) 



where the symbol c?" indicates differentiation of a with respect to e, etc. Note that a and b in equation (A3) are 
Nambu-Keldysh matrices as defined in equation (0). 
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TABLE I. Microscopic material parameters (Y-Ba-Cu-O). 





N f 


Vf 


h/r±\T c 




(K) 


l/(cV cell spin) 


(km/s) 


(meV) 


(mcV) 


92 


4.4 


110 


0.6 


14 



TABLE II. Consistency checks (Y-Ba-Cu-O). 







A,, (0) 
(nm) 


A ± (0) 

w 


dHl 2 /dT\ Tc 
(T/K) 


dH^/dT\ Tc 
(T/K) 


Measured 




- 140 


~ 1 


10 


- -1.8 


s-wave OP 


130 


- 160 


1.1 


-14 


-2.2 


d-wavc OP 


130 


- 230 


7-10 


-55 


-1.9 



TABLE III. Consistency checks (Bi-Sr-Ca-Cu-O). 



Relation: 


Bi-Sr-Ca-Cu-0 


Pb-Bi-Sr-Ca-Cu-0 




Ar-annealed 


02-annealed 


4tt 2 AiA 
he 2 Q± 


- 1.5 


- 5.9 


87re 2 , 


- 1.0 


- 0.20 




- 0.67 


- 0.033 
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FIG. 1. Polar plot of the anisotropic scattering rate, r^ 1 (s,s'), which is normalized here by the 
area for easier visualization. A quasiparticle traveling along the x axis in layer £ (left to right) is 
scattered in a new direction in layer £ ± 1 with a probability proportional to the radial distance on 
this curve. The anisotropy parameter values, 7, shown are: (dashed line), 0.5 (thick line), 1 (dotted 
line), 10 (dash-dotted line). 
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co/A 

FIG. 2. The density of states and the real part of 
the in-plane and c-axis conductivities at T = 0, for a 
scattering rate described by a — = 0.01, and 

several values of the normalized scattering cross section 
a, as a function of frequency. All results are for a d-wave 
pairing interaction unless otherwise noted, (a) Density of 
states normalized by the density of states at the Fermi 
level, (b) The in-plane conductivity normalized by the 
normal state Drude value at u> = 0. The corresponding 
s-wave result (not shown) becomes nonzero at Tiu> /Ao = 2 
and essentially follows the Drude curve, (c) The c-axis 
conductivity normalized by the normal state value. 




co/A o 

FIG. 3. The density of states and the real part of the 
in-plane and c-axis conductivities at T = 0, for a scatter- 
ing rate described by a = Ti/2tttAo = 0.1, and several 
values of the normalized scattering cross section a, as a 
function of frequency. All results are for a d-wave pairing 
interaction unless otherwise noted, (a) Density of states 
normalized by the density of states at the Fermi level, 
(b) The in-plane conductivity normalized by the normal 
state Drude value at u) = 0. (c) The c-axis conductivity 
normalized by the normal state value. 
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FIG. 4. The temperature dependence of the in-plane and c-axis conductivities for a scattering rate 
of a — Tl/2-ktAo — 0.01, in the unitarity limit, as a function of frequency, a) The in-plane conductivity 
normalized by the normal-state Drude conductivity, b) The c-axis conductivity normalized by the 
normal-state value. 
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FIG. 5. The imaginary part of the c-axis conductivity (reactance) at T = for a scattering rate 
described by a = h/2-KTAo — > 0. The s-wave results are denoted by an S. We show 02(01) = Ima ± (u), 
and separately the London term, A/ttuj, and the difference of these two. The constant A was extracted 
from ct ± (uj) by a Kramers-Kronig transformation. Inset: The d- wave result has been rescaled to yield 
the same penetration depth, A±(0), as for the s-wave result (see text). The shaded region indicates 
the range of measured plasma frequencies, u> pc , in cuprate superconductors. 
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